Observables and initial conditions for self-similar ellipsoidal flows 
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Single-particle spectra and two-particle Bose- Einstein correlation functions are determined ana- 
lytically utilizing a self-similar solution of non-relativistic hydrodynamics for ellipsoidally-symmetric, 
expanding fireballs, by assuming that the symmetry axes of the ellipsoids are tilted in the frame of 
the observation. The directed, elliptic and third flows are calculated analytically. The mass depen- 
dences of the slope parameters in the principal directions of the expansion, together with the mass 
and angular dependences of the HBT radius parameters, reflect directly the ellipsoidal properties of 
the flow. 



Introduction — The equations of hydrodynamics re- 
flect the local conservation of matter, momentum and en- 
ergy. These equations are well suited to study problems 
related to flows in various flelds ranging from evolution 
of galaxies in astrophysics to heavy-ion and elementary- 
particle collisions in high-energy physics. The flnding of 
exact self-similar hydro solutions sometimes represents 
essential progress in physics, as the discovery of the Hub- 
ble flow of our Universe or the Bjorken flow of ultrarela- 
tivistic heavy-ion collisions. 

In this Letter we consider the case of a non-relativistic 
hydrodynamical problem with ellipsoidal symmetry. Our 
goal is to demonstrate the influence of initial conditions 
on the final state observables, utilizing an explicit, exact 
and simple analytic solution of fireball hydrodynamics. 
In particular, we attempt to understand the relationship 
between the initial conditions (the ellipsoidal asymmetry 
and the tilt of the major axis) and the final observables. 

The new family of self-similar ellipsoidal solutions — 
A self-similar solution of non-relativistic hydrodynam- 
ics with ideal gas equation of state and a generalized 
(direction-dependent) Hubble flow, a three-dimensional 
ellipsoidal Gaussian density profile and a homogeneous, 
space-indepenent temperature profile has just been found 
in ref. [1]. This solution has many interesting properties, 
e.g. the partial differential equations of hydrodynamics 
are reduced to ordinary differential equations correspond- 
ing to a Hamiltonian motion of a massive particle in a 
non-central repulsive potential. These results correspond 
to the generalization of earlier, data motivated hydrody- 
namical parameterizations and/or solutions of refs. [2-8] 
to ellipsoidal symmetry and non-central heavy-ion colli- 
sions with homogeneous temperature profile. 

Below we generalize the hydrodynamic solution of 
ref. [1] for some fairly wide family of thermodynamically 
consistent equations of state, and calculate analytically 
all the observables of non-central collisions. It allows, in 
principle, to solve an inverse problem, namely, given an 
(in general non-ideal) equation of state, to restore the ini- 
tial conditions from the observables. However, we do not 
aim here to apply directly the new hydro solution to data 
fitting in high-energy heavy- ion physics. In order to reach 
the level of data fitting, generalizations to relativistic flow 



patterns, more realistic equations of state and tempera- 
ture profiles are needed. Some of these generalizations 
seem to be straightforward and are in progress [9,10]. 

Consider the non-relativistic hydrodynamical problem, 
as given by the continuity, Euler and energy equations: 



dtn + V ■ (n-v) = , 
a* v+ (v V)v = -{Vp)/{mn) , 
9t e + V • (ev) = -pV • v , 



(1) 
(2) 
(3) 



where n denotes the particle number density, v stands for 
the non-relativistic (NR) flow velocity fleld, e for the NR 
energy density, p for the pressure and in the following the 
temperature fleld is denoted by T. This set of equations 
are closed by some equation of state (EoS). We have cho- 
sen analytically solvable generalization of NR ideal gas 
equations of state: 



p = uT , e = K(T)nT, 



(4) 



which allow to study the solutions of NR hydro equa- 
tions for any temperature dependent ratio of pressure to 
energy density, p/e = 1/k{T). The EoS (4) are thermo- 
dynamically consistent for any function k(T), as can be 
checked by using the free energy density /(T, n) and the 
relations 



p = n 



,9/ 
dn 



-f: 



^f-T 



dl 
dT 



(5) 



The function k(T) characterizes the p/e ratio for a broad 
variety of materials: e.g. a non-relativistic ideal gas 
yields k{T) — 3/2. Note, that for finite-size systems 
phase transition can occupy certain temperature inter- 
val similar to a crossover. Then one can model such a 
change of the pressure to energy density ratio at phase 
transition from deconfined quark matter to hadronic one 
by means of smooth variations of the values of k{T) in 
certain temperature domain. Note also that it is usual 
to introduce the speed of sound as — dp/de — 1/k(T), 
so we model the change in the equation of state essen- 
tially with the help of a temperature dependent speed of 
sound. 

For reasons of convenience we choose n, v and T as 
the independent hydrodynamic variables. The NR hy- 
dro equations are solved, similarly as it was done for the 



case of NR ideal gas EoS in ref. [1] , by the following self- 
similar, ellipsoidally symmetric density and flow profiles: 
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n(t,r) = no-exp -— - — - — 



v'(i,r') 
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(6) 
(7) 



where the variables are defined in a center of mass frame 
K' , but with the axes pointing to the principal directions 
of the expansion. The time dependent scale parameters 
are denoted by {X,Y,Z) = {X{t),Y{t),Z{t)), the typi- 
cal volume of the expanding system \s V = XYZ, and 
the initial temperature and volume arc Tq = r(fo) and 
Vq = V{to), and no is a constant. The time evolution of 
the radius parameters X, Y, Z and temperature T are 
governed by the ordinary differential equations 



XX = YY = ZZ = 



■d,^^ ^ X Y Z 
^^^^^^^^[x^Y^Z 



m 

= 0. 



(8) 
(9) 



Note that the equation for the time dependence of the 
temperature can be integrated in a straigthforward man- 
ner to find 



Vo 
V 



dT' 

= exp[K(T)-«(To)]exp / ^«(T') , (10) 

J To 



and this equation further simplified, in the case of a tem- 
perature independent ciS 



T = Tn 



l/fv 



(11) 



Observables from the new solution — In order to evalu- 
ate the measurable quantities, any hydrodynamical solu- 
tion has to be supplemented with an additional freeze-out 
criterion, that specifies the end of the hydrodynamical 
evolution. Here we assume sudden particle freeze-out at 
a constant temperature T{tf,r) = Tf where EoS corre- 
sponds approximately to ideal gas {n{Tf) = 3/2). This 
freeze-out condition is reached everywhere at the same 
time in the considered class of exact hydrodynamical so- 
lutions and it is motivated by the simplicity of the results. 
Then, the emission function is proportional to 



S'(t,r',k') oce 



5{t-tf). (12) 



Single particle spectrum — The single-particle spec- 
trum and the two-particle correlation function can be 
evaluated similarly to that of ref. [3,4]: 
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Z=Tf+mXj 
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: Tf + mYf 



Z=Tf+mZj 



(14) 
(15) 
(16) 



Here E = m + k'^/(2m) in the non-relativistic limit we 

are considering, k' = {k'^^k'y^k'^) stands for the momen- 
tum vector in K' , Xf — X{tf), etc. In the spherically 
symmetric case oi X = Y — Z = R, wc recover the 
earlier results [3,4], with (u) = R and T^s = Tf +m{uy. 

The observables are determined in the center of mass 
frame of the collision, K, where the axis points to the 
direction of the beam and the axis to that of the im- 
pact parameter. In this frame, the coordinates and the 
momenta are denoted by x and k. We assume that the 
initial state of the hydrodynamic evolution corresponds 
to a rotated ellipsoid in K. The tilt angle 6 represents the 
rotation of the major (longitudinal) direction of expan- 
sion, from the beam axis r^. Hence the event plane 
is the {r'^,r'^) plane, which is the same as the {rx,rz) 
plane. The (zenithal) angle between directions and 
is 0, while the (azimuthal) angle between the transverse 
momentum k( and the event plane is (f> . 

The ellipsoidal spectrum of eq. (13) generates the fol- 
lowing ^ averaged single-particle spectrum in the K 
frame: 
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where f{v,w) is calculated for \v\ <C 1 and /„(w) = 
i Jq dz cos{nz) exp[w cos(z)] is the modified Bessel func- 
tion of order n (n = 0, 1, ...). For small ellipsoidal asym- 
metries, w -C 1, Io{w) ~ 1 and the effective tempera- 
ture parameter in the transverse direction is the harmonic 
mean of the temperature parameters of the principal di- 
rections of expansion (projected to the transverse plane). 
As T^ >T^ > Ty is expected from the initial conditions, 
we obtain T^ > T^ and Tz<T'^. 

The flow coefficients w„ are defined as 



dh 



dkzktdktdc/) 2'Kdkzktdkt 



1 + 2 Vn cosincj)) 



■ (24) 



Here vi is called the directed flow. V2 the elliptic flow 
and V3 the third flow. The transverse- and longitudinal- 
momentum dependence of the Vn flow components can 
be written in terms of v and w. Assuming that the tilt 
angle 9 or the anisotropy is small, \v\ <C 1, the directed, 
elliptic and third flow components are evaluated as 



V 



1 + 



h{w) 
Io{w) 



Ii{w) v"^ 
2 h{w) 



h{w) h{w) 



Io{w) \lo{w) 



(25) 
(26) 
(27) 



An angular tilt 9 ^ is evidenced by the rise of the 
directed and third flows as a function of rapidity y = 
0.51n[(£' + kz)/{E — kz)] and by a minimum of the el- 
liptic flow at mid-rapidity, see Fig. 1. This and other 
featurc;s are in qualitative agreement with most of the 
data on intermediate- and high-energy heavy-ion colli- 
sions [11-14], suggesting that in non-central collisions the 
dominant longitudinal direction of expansion is slightly 
deviating from the beam direction. A more straightfor- 
ward proof of the ellipsoidal nature of the flow can be 
obtained by determining the mass dependence of the pa- 
rameters , and , cf. eqs. (14-16) and Fig. 2. 

Two-particle correlations The two-particle Bose- 
Einstein correlation function (BECF) is related to a 
Fourier-transform of the emission (or source) function 
S'(i,r',k') of eq. (12), see e.g. refs. [3,4,6]. If the core- 
halo picture [15] is valid, an effective intercept parameter 
A = A(k) = [Nc{]<.)/N{'k)]^ appears, that measures the 
fraction of particles emitted directly from the core. The 
two-particle BECF is diagonal in K' , as 

C{K', q') = 1 + A exp {-q'^^K^ - q'y'R'y' - q'.'R'J') , (28) 
K' = K;2 = 0.5(k;+k'2), (29) 

q' = q;2 = k;-k^ = (g;,g;,g;), (so) 



(31) 
(32) 
(33) 



These radius parameters measure the lengths of ho- 
mogeneity [16]. They are dominated by the short- 
est of the geometrical scales {Xf,Yf,Zf) and the cor- 
responding thermal scales deflned by {Xt,Yt, Zt) = 

"^-(^, ^), generalizing the results of refs. [3,6,7] 

to ellipsoidal flows. The geometrical scales characterize 
the spatial variation of the fugacity term fj,{t, r')/T{t, r'), 
while the thermal scales characterize the spatial varia- 
tions of the Boltzmann term E[^^{t,r')/T{t,r'), both of 
them evaluated at the point of maximal emittivity. In 



the K' frame, cross-terms [17] vanish, -R^^j = 0, if the 
emission is sudden. 

If the particle emission is gradual, but it happens in 
a narrow interval At centered at t / , then the BECF 
can be evaluated using the replacement 6{t — tf) 
(27rAt2)-V2exp[-(t-i/)2/2At2] ineq. (12), ifAtXf « 
Xf , etc. Hence all the previous radius components, in- 
cluding the cross-terms, are extended with an additional 
term SR'^J = 0'if3'jAt'^, where 0' = (k'l + \<i'2)/{E[ + E'^) 
is the velocity of the pair in K' . 

The BECF's are usually given in the side-out- 
longitudinal or Bertsch-Pratt (BP) parameterization. 
The longitudinal direction, riong = n in BP coincides 
with the beam direction. The plane orthogonal to the 
beam is decomposed to a direction parallel to the mean 



transverse momentum of the pair, rout 



and the 



one perpendicular both to this and the beam direction, 
?'sido = ''s • The mean velocity of the particle pair can be 
written in BP as /3 = (/3o,0,/3i), where (3o = Pt ■ Let (p 
denote the angle of the event plane and the mean trans- 
verse momentum of the measured pair. The result is 



C2(K, q) = 1 + A exp - ^ q.q^Rl 

\ s,o,l 

R'^ = R'^ cos^ 4> + Rl sin^ (j) , 
Rl = Rl cos^ 4> + R'^ sin^ cj) + l3^At^, 
Rf = R'J" cos^ 9 + R'^^ sin2 9 + (3fAt\ 



(34) 

(35) 
(36) 
(37) 

Rl^ = [R'^ - i?;2) cos 9 sin 6* cos -I- /?tAA^^ (38) 



{Rl - R'y^) con (t>s\ncj), 
Rl^ = [r!^ - R'^^) cos 9 sin 9 sin < 



where an auxiliary quantity is introduced as 
Rl = R'^ cos^ e + R'^ sin^ 6 . 



(39) 
(40) 

(41) 



These results imply that all the radius parameters oscil- 
late in the K frame. In particular, a <f) dependent os- 
cillation appears in the radius parameters indexed either 
by the side or the out direction, as illustrated in Fig. 
3. These oscillations are similar to those obtained in 
ref. [18], corresponding to ^ = 0. We find that the ra- 
dius parameters indexed by the longitudinal direction de- 
pend also on the zenithal angle 9. A toy model for tilted 
{9 ^ 0) ellipsoidal static pion sources was introduced in 
refs. [19], to understand the (6 dependent oscillations of 
measured HBT radii at AGS energies. In our case, the 
amplitude of the oscillations is reduced for heavier par- 
ticles due to the hydrodynamic expansion, which results 
in a decrease of the lengths of homogeneity with increas- 
ing mass. The oscillations of the radius parameters were 
not related before either to hydrodynamic flow with ellip- 
soidal symmetry and tilt of the major axis or connected 
to the initial conditions of a hydrodynamic expansion. 

A check of the applicability of our hydrodynamic solu- 
tion is that the BECF and the single particle spectrum 



become diagonal (after removing a term of Pi/SjAt'^ from 
all the HBT radius parameters) in the same frame, see 
eqs. (13,14-16) and eqs. (28,31-33). This frame is K' , the 
natural System of Ellipsoidal Expansion or SEE. 

Summary — We have analytically evaluated the single- 
particle spectrum and the two-particle Bose-Einstein cor- 
relation function for a self-similarly expanding, exact, el- 
lipsoidal solution of the non-relativistic hydrodynamical 
equations, assuming a constant freeze-out temperature. 

The parameters of the hydro solution at freeze-out time 
Tf, {Xf,Yf,Zf) and {Xf,Yf,Zf), can be reconstructed 
from the measurement of the single particle spectrum and 
the two-particle correlation functions. The direction of 
the major axis of expansion in the center of mass frame 
of the collision is characterized by the polar angles (9, 0). 
With the exception of all the radius parameters os- 
cillate as function of (p, while the radius parameters , 
, depend also on 0. All the radius parameters de- 
crease with increasing mass, including all the cross terms. 
If 7^ 0, the effective temperature in the transverse direc- 
tion is increased by a contribution from the longitudinal 
expansion. 

The initially more compressed longitudinal and impact 
parameter directions (the and directions) expand 
more dynamically [1,9], that imphes T'^>T'^>T'y. The 
initial time can be identified from the requirement that 
y(io) = 0. The initial conditions for this hydrodynami- 
cal system can he uniquely reconstructed from final state 
measurements. The function k(T) in the EoS influences 
only the time evolution of the scales (X, F, Z) and tem- 
perature T. So, for a given EoS, one can uniquely recon- 
struct the initial conditions of hydrodynamic evolution 
from final ones. 

We have deliberately chosen the presentation as sim- 
ple as possible, which limits the direct applicability of our 
results in high-energy heavy-ion collisions only to suffi- 
ciently small transverse momentum, pt <^ m, at mid- 
rapidity. But the scheme permits generalization in many 
points and still the qualitative features of our results may 
survive even in the relativistic regime. Generalizations 
to some relativistic flows have been described in ref. [10]. 
Ref. [9] includes another extension to an arbitrary, inho- 
mogeneous, ellipsoidally symmetric initial temperature 
profile, which does not change the time evolution of the 
scale parameters. It turns out that each of these gener- 
alizations is essentially straightforward. 

Although we considered a non-relativisitic problem, 
our results provide generic insight into the time evolution 
of non-central heavy-ion collisions and relate initial con- 
ditions to final-state observables in a simple and straight- 
forward manner impossible before. 
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FIG. 1. The directed, elliptic and third flows vi, V2, 
V3 are illustrated, respectively with solid, long-dashed and 
short-dashed lines, as a function of rapidity, for m — 940 
MeV, = 200 MeV, Ty = 150 MeV, = 700 MeV, at a 
fixed kt = 500 MeV and 9 = 7r/5, see eqs. (25-27). 
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FIG. 2. The linear mass dependence of the effective tem- 
peratures in the transverse directions, T^, Ty and their (har- 
monic) average T^s for non-central heavy-ion collisions, if 
Tf = 120 MeV, X = 0.4, and Y = 0.25 and 8 = 0. 
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FIG. 3. The out-side and the side-long cross terms are plot- 
ted as a function of the polar angle 0, for R'^ = 5 fm, ii^ = 4 
fm, 7?2 = 8 fm and 9 = 7r/5. Note that R^i oc sm{4>) while 
i?os oc sin(2(^). See eqs. (39) and (40). 



